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Abstract An analytical approximation of a nonlinear force-free magnetic field 
(NLFFF) solution was developed in Paper I, while a numerical code that per- 
forms fast forward-fitting of this NLFFF approximation to a linc-of-sight mag- 
netogram and coronal 3D loops has been described and tested in Paper II. Here 
we calculate the free magnetic energy Ef roo = £"n — £p, i.e., the difference of the 
magnetic energies between the nonpotential field and the potential field. A sec- 
ond method to estimate the free energy is obtained from the mean misalignment 
angle change Afj, = (J,p — /Un between the potential and nonpotential field, which 
scales as E{ ree /Ep tan 2 (A/i). For four active regions observed with STEREO 
in 2007 we find free energies in the range of qw = {Ef ICC /Ep) » 1% — 10%, 
with an uncertainty of less than ±2% between the two methods, while the free 
energies obtained from 11 other NLFFF codes exhibit a larger scatter of order 
»s ±10%. We find also a correlation between the free magnetic energy and the 
GOES flux of the largest flare that occurred during the observing period, which 
can be quantified by an exponential relationship, -Fgoes oc ex P (<7free/0.015), 
implying an exponentiation of the dissipated currents. 

Keywords: Sun: Corona — Sun: Magnetic Fields 



1. Introduction 

The free magnetic energy is the maximum amount of energy that can be released 
in an active region of the solar corona, such as during a solar flare, a filament 
eruption, or a coronal mass ejection (CME). Thus, it is important to design 
reliable methods and tools that can calculate the amount of free energy in 
order to quantify the energy budget in a catastrophic energy release event, as 
well as for estimating upper limits in forecasting individual events in real-time. 
Traditionally, the free energy is calculated by computing the nonpotential field 
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Bn(x) with a numerical nonlinear force- free field (NLFFF) code and a potential 
field Bp(x) for the same photospheric boundary data B(x, y, z p hot), so that 
the difference of the magnetic field energy density integrated over a volume V 
encompassing the active region of interest can be quantified as 



This standard method may not necessarily reflect the correct amount of max- 
imum free energy released during a solar flare, since the magnetic field in the 
photospheric boundary B(;r, y, %>hot) may change during a flare (e.g., see mea- 
surements by Wang et ai, 1994, 2002, 2004, 201 3; Wang, 1997, 2006; Wang and 
Liu, 2010). Another problem with NLFFF codes using the photospheric vector 
field is the non- force- freeness of the lower chromosphere (Metcalf et al, 1995; 
DeRosa et al., 2009), which however, can be ameliorated by preprocessing the 
magnetic boundary data, using chromospheric field measurements (e.g., Metcalf 
et ai, 2005; Jing et ai, 2010), or by a multigrid optimization that minimizes a 
joint measure of the normalized Lorentz force and the divergence of the magnetic 
field, as proposed by Wiegelmann (2004) and applied by Jing et ai, (2009). A 
quantitative comparison of NLFFF computation methods, however, revealed a 
substantial scatter of free energies in the order of « ±10% (of the potential 
energy), depending on the numeric code, the boundary specifications, and the 
spatial resolution (Schrijver et ai, 2006; DeRosa et al., 2009). Alternatively, 
some studies show that the free energy is better estimated by the minimum- 
energy state above the linear force-free field with the same magnetic helicity 
(Woltjer, 1958; Regnier and Priest, 2007). Since NLFFF calculations are rather 
computing-intensive for forward- fitting tasks to coronal constraints, which re- 
quires many iterations (for an overview and discussion of different numerical 
methods see recent reviews by Aschwanden (2004) or Wiegelmann and Sakurai 
(2012)) or faster non-numerical methods are desirable. Some proxy of the active 
region's free magnetic energy has been defined based on the twist and magnetic 
field orientation near the neutral line (e.g., Falconer et ai, 2006, 2011). 

Free magnetic energies have been calculated for a variety of solar phenomena, 
for instance for the evolution of the free magnetic energy during flux emergence 
and cancellation, using NLFFF codes (e.g., Fang et ai, 2012), for the evolu- 
tion of active regions (e.g., Kusano et ai, 2002), for helmet-shaped streamer 
configurations (Choe and Cheng, 2002), or for breakout CMEs, using MHD 
simulations (e.g., DeVore and Antiochos, 2005). Theoretical studies quantify 
the evolution of free magnetic energy for dipolar (Regnier, 2009) and quadrupo- 
lar magnetic configurations with a null-point, using different force- free models 
(Regnier, 2012). 

An analytical approximation of a nonlinear force- free magnetic field (NLFFF) 
solution was developed in Paper I (Aschwanden, 2012), while a numerical code 
that performs fast forward-fitting of this NLFFF approximation to coronal 3D 
loops and a linc-of-sight magnctogram has been described and tested in Pa- 
per II (Aschwanden and Malanushenko, 2012). In this Paper III of the series 
we are concerned with the calculation of the free energy, generally defined by 
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the difference of the magnetic energies between the nonpotential field and the 
potential field (Equation (1)). We calculate free energies from simulated data 
(from Paper II). from the analytical NLFFF solution of Low and Lou (1990), 
from active region NOAA 10930 during an X3.4 flare modeled by Schrijver 
et al., (2008) and Malanushenko et al, (2012), and from stereoscopically tri- 
angulated loops observed with STEREO (Aschwanden et al., 2012). Section 
2 describes the analytical treatment, Section 3 the application to simulated 
datasets and observations, Section 4 contains a discussion, and Section 5 the 
conclusions. 



2. Analytical Formulation 
2.1. The NLFFF Approximation 

In Paper I (Aschwanden, 2012) we derived an analytical approximation of a non- 
linear force-free field (NLFFF) solution, which fulfills Maxwell's divergence-free 
equation (V • B — 0) and the force-free equation (V x B) = a(x)B with second- 
order accuracy (of the force- free parameter a). The analytical approximation 
can be specified by a radial field B r and an azimuthal field component B v , 

BJr, 9) = Bj ( 4^1 ~ o— , (2) 

V ' J \r 2 J (l + b 2 r 2 sm 2 9) ' V ' 



/ f\ \ ( d 2 \ br sin 9 , . 

B v (r, 6) = Bj i — \ j—r^-. 2 , 3) 
\r 2 J (1 + b 2 r 2 sin 9) 

B e (r,9)^0, (4) 

2bcos6> 

tt ( r > e > ~ 77TT2~2 ■ 2 fl x » ( 5 ) 

(1 + o^H sin 9) 

where (r, (p, 9) are the spherical coordinates of a single magnetic field component 
(Bj,Xj,yj, Zj,bj) with a unipolar magnetic charge Bj that is buried at position 
(xj,yj,Zj), has a depth d = 1 — [x 2 — y 2 — Zj} 1 ^ 2 , a vertical twist a = 2b j, and 
r = [(x — Xj) 2 + (y — yj) 2 + (z — Zj) 2 ] 1 / 2 is the distance of an arbitrary coronal 
position (x,y, z) to the subphotospheric location (xj,yj,zj) of the buried mag- 
netic charge. The force-free parameter a is expressed in terms of the parameter 
b, which quantifies the number A" tw ist of full twist turns over a (loop) length L, 

27rAt w i s t 

b= — — (6) 

This analytical approximation is divergence-free and force-free to second-order 
accuracy in the parameter (b r sin 9), which is approximately proportional to the 
force- free parameter a as defined by Equation (5). 
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A general magnetic field configuration can be composed by a superposition 
of N m twisted magnetic field components, 

»„ 

B N (x)=^B J (x) , (7) 

i=i 

which also fulfils the divergence-free and force-free condition with second-order 
accuracy in a (or b). 

2.2. The Free Magnetic Energy of a Single Twisted Component 

Let us calculate now the free magnetic energy cLEf roo (x) at location (x) for the 
field resulting from a single twisted (buried) magnetic charge as defined by Equa- 
tions (1) to (6). Since the radial B r and azimuthal components B v are always 
orthogonal to each other (Figure 1), we can calculate the total nonpotential 
magnetic field strength £?n at every given point (r, if) simply from the sum of 
the squared components B r and B v , 



B N = (Bt + B<)< (8) 



^(1 + b 2 r 2 sin 2 i 



while the field strength Bp of a potential field corresponds to the radial compo- 
nent B r (of a single buried magnetic charge, Equation (2)), 

cP\ 1_ 

r 2 J (1 + b 2 r 2 sin" 



Bp = Br = S . )_____, (9) 



and thus the free energy dE{ lee (x) is just the magnetic energy associated with 
the azimuthal field component B^, with Equations (8) and (9), 

dE hcc {x) = dE N {x)-dE P {x) = — [S r (x) 2 + J B v (x) 2 ]-^- J B r (x) 2 = -^B^(x) 2 . 

(10) 

This definition of the free magnetic energy dEf iee (x) fulfills the following condi- 
tions: 

1. Positivity constraint: The free energy is positive for every nonpotential field 
at any location x, i.e., dEf rco (x) > 0, since B 2 (x) > and dEf rcc (x) oc B 2 (x) 
according to Equation (10). 

2. Additivity of energies and orthogonality of magnetic field components: The 
nonpotential energy corresponds to the sum of the potential energy and the 
free energy dE^(x) = dEp(x) + dEf rcc (x) according to Equation (1) and 
Figure 1. Since the energies scale with the square of the magnetic field compo- 
nents, i.e., dE^(x) oc Bn(x) 2 , dEp(x) oc B r (x) 2 , and cLEf roo (x) oc B v (x) 2 , the 
Pythagoras theorem can be applied and it follows that B ri B v , and B^ form a 
rectangular triangle, with B v perpendicular to the potential field component 
Bp = B r . 
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(x,y.z) 

Figure 1. Diagram of the nonpotential field vector Bfj, which is composed of the two orthog- 
onal components of the azimuthal field vector B v and the radial field vector B r , subtending 
an angle 8. 

3. Positive scaling with force-free parameter: The free energy quantifies an ex- 
cess of nonpotential magnetic field energy compared with the potential field 
energy, and thus should have a positive scaling with the force-free parameter 
a. Since the free energy scales proportional to the square of the azimuthal 
magnetic field component (Equation (fO)), we have a positive scaling, which 
is approximately t£Ef rcc (x) oc b 2 oc a 2 (Equations (3) and (5)). 

4. Potential field limit: The free energy vanishes asymptotically (c£Ef rcc (x) i-> 0) 
with vanishing force-free parameter \a\ > or \b\ i— > 0, with the potential 
field being the limit, dE^(x) H> dEp(x). 

5. Finiteness of nonpotential energy: The total nonpotential magnetic energy 
integrated over an arbitrary large height r converges to a finite value, E(r t— > 
oo) = E max . We can prove the finiteness of the potential energy for a single 
magnetic charge, which has a square-dependence of the magnetic field, B(r) oc 
r~ 2 , yielding a 4th-power dependence of the magnetic energy dEp(r) oc 
B{r) 2 oc r~ 4 , and thus a 3th-power dependence for the integrated magnetic 
energy, -Ep^ot °c J dEp(r)dr oc r~ 3 . For a finite amount of twist and a finite 
number of magnetic sources, it can be shown that the integral of the resulting 
nonpotential energy is also finite. 

2.3. The Free Magnetic Energy of Multiple Twisted Components 

While the foregoing definition of the free energy is calculated for a single twisted 
(buried) magnetic charge, how can it be generalized for a superposition of an 
arbitrary number of magnetic charges, as defined in Equation (7)? The sum of 
the magnetic field contributions from each buried magnetic charge component 
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Bj (x) add up to the nonpotential field vector Bn (x) , which can be decomposed 
into two orthogonal components B||(x) and Bj_(x) in every point of space (x), 



B N (x) = ^B J -(x)=B||(x)+B ± (x) , (11) 

where the parallel component is aligned with the potential field direction, B|| || 
Bp, and the perpendicular component is orthogonal to the potential field di- 
rection, B^ _L Bp. For a single (twisted) magnetic component the parallel 
component B|| is identical with the radial component B r , and the perpendicular 
component Bj^ is identical with the azimuthal component B^, (Figure 1). The 
three magnetic field components Bu = Bp, B±, and £?n are then associated each 
with one of the three energy components, 

d£kee(x) = (1/87T)B2 ( x ) 

dEp(x) = (l/87r)Sj|(x) (12) 
dE^ (x) =(l/87r)[Sf(x) + Bi(x)] . 

Alternatively, the magnetic energy dEp (x) of the potential field can be computed 
by using current- free magnetic field components (Bj,Xj,yj,Zj,bj = 0) straight- 
forward with Equation (1), and the magnetic energy dE-^fjx.) of the non-potential 
field with the current components (Bj, Xj, yj, Zj, bj ^ 0) with Equation (1) also, 
which yields the free energies (Ef Tee = E^ — Ep), after volume integration. Both 
methods are fitting the same line-of-sight component of the photospheric bound- 
ary B z {x, y, Zphot) given by the magnetogram, while the non-potential magnetic 
field affects the transverse field components B x {x, y, z p hot) and B y (x, y, z p hot) 
that are not used as a boundary condition in our forward-fitting method. 

2.4. Free Energy Estimated from Misalignment Angle with Loops 

An alternative method is to estimate the free energy in an active region from 
the misalignment angle /x between the potential field and the observed coronal 
loops, or the best-fit nonpotential field. Using the coronal loops as a proxy for 
the nonpotential field, this would provide a very fast method that requires only 
the computation of a potential field, supposed we have the 3D coordinates of 
coronal loops, e.g., from stereoscopic triangulation. Since the mean azimuthal 
field component is B± = £?n sin (A/x) and the mean radial (potential) field 
component is B\\ = -Bn cos (A/x), the nonpotential energy ratio = E-^/Ep 
follows then directly, using the definition of Equation (11), for the free energy, 

T? B 2 + B 2 
q, = ^= ± B , " = 1 + tan 2 (Ait) , (13) 

where the relative misalignment angle A/x is defined as the difference between 
the (median) potential /xp and (median) nonpotential field directions /in, 

Afx = /x P - /x N . (14) 
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Figure 2. Forward-fitting of simulated nonpotcntial field data (cases N7-N9). Each panel 
shows the linc-of-sight magnetogram (grey ) , simulated coronal loops (the targets of the forward 
fit) (blue), and magnetic field lines of a theoretical model (red), either the potential field 
constrained by the line-of-sight magnetogram (left panels), or the forward fit of the NLFFF 
approximation (right panels). 
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Figure 3. Forward-fitting of cases N10-N12. Representation similar to Figure 2. 



Ideally, if the NLFFF forward-fitting code matches the coronal loops perfectly 
(with /l«n ~ 0), and thus the relative misalignment angle of the potential field to 
the loops can be used, A/z ~ //p. However, in reality there is always a significant 
difference between the best-fit NLFFF solution and the observed loop data, 
either due to stereoscopic measurement errors or due to an additional field 
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misalignment that cannot be described with the particular parameterization of 
our NLFFF approximation. 

3. Numerical Tests and Results 

We perform now tests of the calculation of the free magnetic energy using the 
definitions given in Equations (12) and (13) for four different datascts, using 
(1) simulated data produced by our analytical NLFFF approximation, (2) the 
NLFFF solution of Low and Lou (1990) with a known exact analytical solution, 
(3) active region NOAA 10930 during an X3.4 flare modeled by Schrijver et 
at, (2008) and Malanushenko et ai, (2012), and (4) stereoscopic data from four 
observed active regions. 

3.1. Tests with Simulated Data (P1-N12) 

We simulated six cases of potential fields (shown in Figure 3 of Paper II) , and six 
cases with nonpotential fields (Figures 2 and 3), which are similar to the cases N7 
to N12 in Paper II, except that we reduced the amount of twist by a factor of five, 
in order to make them more comparable with observations of real active regions, 
which have free energies of order < 10%. We are using the parameterization of 
our analytical NLFFF approximation, described in Paper I, containing 1, 2, 4, 
and 10 magnetic source components with variable vertical twist, labeled as cases 
P1-P6 and N7-N12. We integrate the nonpotential magnetic energy E^ in a box 
that covers the displayed field-of-view centered at the center of the solar disk 
and has a height range of h = 0.15 solar radii above the photosphere. The values 
of the potential energies Ep and the ratio of the nonpotential to the potential 
energy of the model (g m odci), or of the fit (gg t = E-^/Ep), are listed in Table 1. 
The free energy is Ef Tee = £n — Ep = Ep(q^ — 1). 

Table 1 demonstrates several results. First, the misalignment angle between 
the forward- fitted nonpotential field and the simulated loops is in all 12 cases 
(P1,...,N12) significantly smaller (/L*n = 1.4°±0.8°) than the initial potential field 
misalignment angle fip = 9.6°±8.5°, which confirms a satisfactory convergence of 
the forward-fit to the simulated target loops. Second, all ratios of nonpotential 
to potential energies E-^/Ep are larger than one, which means that the free 
energy E{ rcc = E^ — Ep is always positive with our definition given in Equation 
(12). Third, all 12 simulated cases agree in the nonpotential energy ratio with 
the simulated input model, i.e., qat = q m ode\ = E^/Ep with an accuracy of 
less than 10~ 5 , which also confirms the perfect convergence of the forward- fit 
algorithm. The nonpotential-field cases (N7-N12) have free energies in the range 
of q&t = E^/Ep = 1.010—1.163, or up to 16% of the potential-field energy, which 
are also retrieved with an accuracy of better than 10 -5 . We compare also the 
nonpotential energy ratios calculated from the forward-fit (qat) an d estimated 
from the median misalignment angles (q^; Equation (13)) and find an agreement 
of q^/qfit = 1-02 ± 0.05 between the two methods. 

We test also the finiteness of the nonpotential energy. In Figure 4 we plot 
the height-integrated total nonpotential energies E^(r) as a function of the 
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Table 1. Magnetic energy calculations arc listed for 20 cases, including six poten- 
tial-field cases (P1-P6), six nonpotential-field cases (N7-N12), the Low and Lou 
(1990) case (LI, L2), the Schrijver et al., (2008) case (SI, S2), and four stereoscop- 
ically observed active regions (A, B, C, D), specified by the number of magnetic 
charges (7V m ), the fraction of magnetic energy captured by the model (qE), the 
number of simulated loops TV; , the potential energy Ep (corrected by the factor 
<?£;), the median misalignment angle fip of the potential field, the median misalign- 
ment angle fi^ after forward-fitting of the NLFFF model, the predicted energy 
ratio q M = (l.+tan 2 (Aft)) based on the misalignment angle change Aft = ftp — ftjv, 
the forward-fitted nonpotential energy ratio q^t = En/Ep, and the volume-inte- 
grated nonpotential magnetic energy ratio of the model q mo dcl = ^modol/^P; with 
values computed by Anna Malanushenko (private communication, 2012) using the 
Low and Lou (1990) data ( a ), the value of 3D-fits labeled as II. b Tables 3 and 
4 of Malanushenko et al, (2012) ( b ), and the value of the Wh + code with the 
smallest misalignment angle ft = 24° in Table 1 of DeRosa et al., (2009) ( c ). 



Case N m (q E ) A loop ftp ft N E P q^ <j fit ijmodol 

10 32 
(erg) 



PI 


1(1.000) 


49 


0.5 


0.0 


4.65 


1.000 


1.000 


1.000 


P2 


2(1.000) 


121 


1.8 


1.2 


5.84 


1.000 


1.000 


1.000 


P3 


4(1.000) 


121 


2.1 


1.4 


4.27 


1.000 


1.000 


1.000 


P4 


10(1.000) 


121 


2.5 


1.3 


14.6 


1.000 


1.000 


1.000 


P5 


10(1.001) 


121 


3.7 


1.4 


12.9 


1.002 


1.000 


1.000 


P6 


10(0.998) 


121 


3.3 


1.9 


6.80 


1.001 


1.000 


1.000 


N7 


1(1.000) 


49 


18.8 


0.2 


4.65 


1.114 


1.010 


1.010 


N8 


2(1.000) 


121 


10.4 


1.2 


5.84 


1.026 


1.009 


1.009 


N9 


4(1.000) 


121 


13.0 


1.8 


4.27 


1.039 


1.016 


1.016 


N10 


10(1.000) 


121 


17.4 


3.1 


14.6 


1.065 


1.082 


1.082 


Nil 


10(1.001) 


121 


17.2 


1.1 


12.9 


1.084 


1.113 


1.113 


N12 


10(0.998) 


121 


25.3 


1.9 


6.80 


1.188 


1.163 


1.163 



LI 100(0.913) 133 14.7 4.8 0.246 1.030 1.023 1.129" 

L2 100(0.910) 35 3.9 1.6 0.028 1.002 1.001 1.091" 

SI: 2006/12/12 200(0.860) 331 37.3 14.4 18.3 1.179 1.112 1.21 ±0.05 6 

S2: 2006/12/13 200(0.878) 98 27.3 13.4 17.8 1.061 1.104 1.08 ± 0.01 6 



A: 2007/04/30 


100(0.854) 


200 


23.8 


21.4 


12.7 


1.002 


1.006 


1.030 c 


B: 2007/05/09 


100(0.807) 


70 


21.2 


17.9 


1.08 


1.003 


1.023 




C: 2007/05/19 


100(0.784) 


100 


41.3 


22.1 


3.47 


1.121 


1.085 




D: 2007/12/11 


100(0.810) 


87 


22.1 


14.7 


7.83 


1.017 


1.044 





height limit r for the six cases N7-N12 (Table 1) and find that each one follows 
approximately an exponential height dependence (dashed curves in Figure 4), 

£4.(r)w£U[l-fflp(-^)], (15) 

which has the finite limit E(r i— > oc) = E max , and thus confirms the convergence 
of the code. In our calculations we generally use a height limit of h max = 0.15 
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0^ , , , , i , , L , i , , , , 

1.00 1.05 1.10 1.15 

Altitude r[solar radii] 



Figure 4. Height dependence of the nonpotcntial energy Ej^(h) integrated over a volume from 
height r = 1.0 to r = r max with r max = 1.15 solar radii for the six-nonpotcntial cases (N7-N12 
in Table 1). An exponential function is fitted in the upper half height range r > 1.075 (dotted 
line) . 



solar radii, which corresponds to about two density scale heights of T e = 1.0 MK 
plasma, or four emission measure scale heights. 

3.2. Tests of the Low and Lou (1990) NLFFF Solution 

An analytically exact solution of a NLFFF model was derived by Low and 
Lou (1990), described also in Malanushenko, Longcope, and McKenzie (2009). 
The particular solutions we are using are defined by the parameters (a=0.6, 
n=2.0) in case LI, and (a=0.01, n=1.0) in the case L2, where a is the Grad- 
Shafranov constant, n is the harmonic number of the Legendre polynomial, and 
additional parameters are the depth I of the source below the photosphere and 
the inclination angle of the axis of symmetry). Forward- fits of our analytical 
NLFFF approximation to the exact NLFFF solution of Low and Lou (1990) 
are shown in Figure 5 (top and middle), where we used a computation box 
of (100 x 100 x 75), with a pixel size of Ax = 0.002 solar radii and a height 
range of h max = 75 Ax « 0.15 solar radii. The forward-fit was accomplished by 
using the line-of-sight magnetogram at a planar surface {B z {x,y,z = 1) and 
a set of Ni = 133 (case LI) and -/Vi oop = 35 (case L2) target field lines that 
mimic coronal loops. The misalignment angle is reduced from /xp = 14.7° to 
/^n = 4.8° for case LI, and from /ip = 3.9° to /in = 1-6° for case L2, so the 
forward-fitting reduces the misalignment by about a factor of 3. The resulting 
nonpotcntial energy ratios are listed in Table 1, yielding free energies of 2.3% for 
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case LI, and 0.1% for case L2, respectively. These free energies are significantly 
below the theoretical values calculated in Malanushenko et ai, (2009), where 
values of 13% and 9% are quoted (listed under g mo doi in Table 1), although the 
misalignment of the forward-fitted field is quite satisfactory. The reasons for 
the mismatch in the free energy for this particular case is not fully understood, 
since the convergence behavior of our code seems to be no problem in case of 
unique solutions (such as the simulated cases N7-N12). We suspect that the 
parameterization of our NLFFF approximation, which consists of a number of 
buried point-like magnetic sources, is not adequate or suitable to represent the 
analytical Low and Low (1990) solution, which consists of extended, smooth 
magnetic distributions with elliptical shapes, parameterized in terms of Legendrc 
polynomials. It is conceivable that the representation of Legendre polynomials 
by (spherically symmetric) point sources leads to clustering of closely-spaced 
point sources with cancelling of the nonpotcntial (azimuthally twisted) field 
components. 

3.3. Tests of the case of Schrijver et at, (2008) 

The magnetic field of active region NOAA 10930, observed with TRACE, Hin- 
orfe/XRT, and Hinode/SOT on 2006 December 12, 20:30 UT before a GOES- 
class X3.4 flare (case SI), and on 2006 December 13, 03:40 UT after the flare (case 
S2), has been extensively modeled with NLFFF codes (Schrijver et ai, 2008; 
Malanushenko et ai, 2012). We show a forward-fit of our NLFFF approxima- 
tion in Figure 5 (bottom panels) to a set of loops (i.e., closed field lines that 
are randomly chosen from a NLFFF solution computed by Malanushenko et 
ai, 2012). The accuracy of the forward- fitting depends mostly on the number 
of magnetic field components N m , but generally reaches asymptotically a flat 
plateau for N m > 100 (Aschwanden et al, 2012; Figure 10). For this magnetically 
very complex active region we needed N m = 200 magnetic source components 
to reach the plateau, while N m < 100 was sufficient for all other cases. Here, the 
median misalignment angle of fip = 37.3° for the potential field was reduced by 
about a factor of 2.6 to /in = 14.4° for the best-fit nonpotcntial NLFFF solution 
(case SI in Table 1). We measure a potential-field energy of Ep = 1.83 x 10 33 
erg before the flare (SI), and E P = 1.78 x 10 33 erg after the flare (S2), so a small 
difference of « 2.6%. For the nonpotential magnetic energy ratio we measure 
qat = E^/Ep = 1.112 before the flare and gg t = 1-104 after the flare, and 
similar values with the misalignment method, i.e., = 1.179 and q^ = 1.061. 
Thus the total nonpotential energy decreases by AE^ = 0.7 x 10 32 erg (w 5% 
according to the forward-fit method), or by AE^ = 2.7 x 10 32 erg (w 13% 
(according to the misalignment method). 

From the same observing times the free energy ratio was measured with 14 
different NLFFF codes in Schrijver et ai, (2008; Table 1 therein), which yield 
energy ratios of E^/Ep : f rcc = 1.05 ± 0.05 before the flare, and E^/Epj ICC = 
1.16 ± 0.14 after the flare, if we average all methods with equal weight. How- 
ever, the most reliable NLFFF method among them, according to a quality 
assessment by visual inspection of five magnetic features seems to be the Whp p 
NLFFF code, which yields an energy ratio of E^/Ep tPYe = 1.32 before the flare 
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and E-^/Ep, plc = 1.19 after the flare, so a decrease of 13% in the free energy, 
corresponding to a drop of AEf ree w 3 x 10 32 erg in free energy, similar to our 
measurement with the misalignment method (£f rcc = 2.7 x 10 32 erg or 13%). 
Similarly, Malanushenko et ai, (2012; cases with 3D fits labeled as II. b in Table 
1 therein) calculates energy ratios of E-^/Ep = 1.21 ± 0.05 before the flare and 
En/Ep = 1.08 ± 0.01 after the flare, corresponding to a drop of 13% in the free 
energy. 

Thus all three studies agree with a drop in free energy, by an amount of 
5% — 13% according to our two methods, versus « 13% for the most reliable 
NLFFF codes, while lower free energy values and larger misalignment angles 
result for the other NLFFF codes (Table 1 in Schrijvcr et al, 2008). Thus our 
two methods appear to be commensurable with the most reliable NLFFF codes. 

3.4. Tests with Stereoscopic Observations 

We calculate now the free magnetic energy for four active regions (Table 1) that 
have all been observed with STEREO and were subjected to previous magnetic 
modeling, which we label as active regions A, B, C, and D. All four active regions 
have been subject of potential-field modeling using stereoscopic data, including 
potential-field stretching (Sandman et al, 2009), buried magnetic charges (As- 
chwanden and Sandman, 2010), or buried dipoles (Sandman and Aschwanden, 
2011), and nonlinear force-free field modeling (Aschwanden et al, 2012). 

Active region NOAA 10953, 2007 April 30 (A) has been observed during a 
flare energy build-up phase and a flare trigger by magnetic reconnection at 
a 3-D null point of a separatrix surface was identified (Su et al, 2009). The 
magnetic modeling of this active region has been scrutinized with 11 different 
NLFFF codes, using SoHO/MDI, ffinode/SOT-SP, XRT, and STEREO/EUVI 
data (DeRosa et al, 2009; Su et al, 2009). For this active region (shown in 
Figure 6, top panels) we measure a potential field energy of Ep = 1.27 x 10 33 
erg on 2007 April 30, while Su et al, (2009) measure E P = 0.96 x 10 33 erg on 
2007 May 2, which agrees within « 30%, over a time difference of two days. For 
the nonpotential energy ratio we find E-^/Ep = 1.006 on 2007 April 30, while Su 
et al., (2009) find w 1.1 on 2007 May 2, in the central core of the active region, a 
few hours before a GOES-class B3.8 and C8.5 flare. Their higher value could thus 
be attributed to flaring activity. Extensive NLFFF modeling was carried out for 
2007 April 30 at 22:24 UT, using Hinode vector magnctograph data (DeRosa et 
al., 2009). The free energies obtained from 11 different NLFFF codes scatter in 
the range of E N /E P = 0.87 - 1.24, with a value of 1.03 for the Wh + code with 
the smallest misalignment angle of (fi = 24°) with respect to the stereoscopically 
triangulated loops. Including the uncertainties of the boundary conditions in the 
NLFFF code, a self-consistent NLFFF solution with a nonpotential energy ratio 
of E^/Ep « 1.08 was obtained, with a potential-field energy of Ep = 0.84 x 10 33 
erg (Wheatland and Leka, 2011). 

Part of the discrepancy could possibly be explained by the different ficld-of- 
view, since the Hinode field-of-view used in DeRosa et al, (2009) covers only 
the central part of the active region, while our field-of-view encompasses the 
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entire active region. Thus we calculated the magnetic energy in the same Hinode 
ficld-of-view as used in DeRosa et ai, (2009) but find almost identical energy 
ratios. Our lower value than obtained with the other NLFFF codes could also be 
attributed to an underestimation of the twist (and thus nonpotcntial energy) in 
the core of the active region, where stereoscopic loop triangulation is very sparse 
due to confusion of loops with "moss" background. If this is true, we generally 
expect that the avoidance of twisted core structures leads to a stereoscopic 
undersampling bias, resulting into lower estimates of the free magnetic energy. 
We also have to keep in mind that the free energy ratio for this active region 
is the lowest among the four active regions, and thus has the largest relative 
uncertainty 

Active region NOAA 10953, 2007 May 9 (B) was subject to the first 3D re- 
construction with STEREO (Aschwanden et ai, 2008a), stereoscopic electron 
density and temperature measurements (Aschwanden et ai, 2008b), and in- 
stant stereoscopic tomography and DEM modeling (Aschwanden et ai, 2009). 
This active region exhibits the simplest bipolar magnetic configuration among 
all four analyzed active regions and we find a nonpotcntial energy ratio of 
E N /E P = 1.023 (Table 3), i.e., a free energy ratio of w 2.3%. 

Active region NOAA 10953, 2007 May 19 (C) has exhibited multiple filament 
eruptions in the complex and highly nonpotential magnetic configuration during 
2007 May 19 (Li et al, 2008; Licwcr et al, 2009; Hara et al, 2011). Some 22 
GOES B-class and 2 GOES C-class flares were detected during the observing 
period (Li et al, 2008). A filament eruption, accompanied by a B9.5 flare, coronal 
dimming, and and EUV wave was observed and traced with 3D reconstruction 
after 2007 May 19, 13:00 UT (Licwcr et al, 2009). The associated EUV dimming 
and EUV wave caused by the filament eruption was also analyzed (Attrill et 
al, 2009). Plasma motion and heating up to T c — 9 MK was observed for the 
same flare around 13:00 UT (Hara et al., 2011). For this active region, which 
we analyzed at 12:47 UT shortly before the flare (shown in Figure 6, middle 
panels), we found the largest amount of free energy (9%), i.e., a nonpotcntial 
energy ratio of E-^/Ep = 1.085 (Table 1), which clearly is associated with the 
filament eruption and flaring activity after 13:00 UT. 

Active region NOAA 10978, 2007 December 11 (D) appears also to have a 
dominant bipolar structure (shown in Figure 6, bottom panels), but some ap- 
parent currents along the central neutral line have been modeled with a flux- 
insertion method (Alex Engcll and van Ballegooijen; private communication 
2012). For this active region we found a moderate amount of free energy (4%), 
i.e., a nonpotential energy ratio of E^/Ep = 1.044 (Table 1), which is likely to 
be stored in the flux rope or filament above the central neutral line. 

3.5. Correlation of Free Energies with Flare Fluxes 

It was noted in earlier work that the level of flaring activity in these four active 
regions is positively correlated with the nonpotentiality of the magnetic field, 
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Figure 7. GOES soft X-ray light curves of the 0.5-4 A (upper curve) and 1-8 A channel 
(lower curve) during the time of stereoscopic triangulation and magnetic modeling of the four 
analyzed active regions. The flare peak and prcflarc background levels are indicated with a 
step function (thick solid lines). The four active regions are in order of increasing GOES flux 
(adapted from Aschwanden and Sandman, 2010). 

evaluated by the average misalignment angle between a magnetic potential-field 
model and the observed 3D loop coordinates. This correlation was interpreted 
in terms of a relationship between electric currents and plasma heating (As- 
chwanden and Sandman, 2010). Since the misalignment angle is a measure of 
the magnetic nonpotentiality, we expect that there should also be a correlation 
between the flaring activity level and the free magnetic energy in a flaring active 
region. The GOES 0.5-4 and 1-8 A light curves are shown in Figure 7. We 
plot the (preflare background-subtracted) GOES 0.5-4 A fluxes /goes of the 
largest flare that occurred during the observing period of an active region versus 
the free energy ratio qB,hee = Es.free/ Ep in Figure 8 (top panel). The preflare 
background-subtracted GOES fluxes Fgoes show a correlation that can be fitted 
with an exponential function, 



( /goes\ 


( <7freeA 


V fo J 


« exp 


V 9o / 



(16) 



with the constants go = 0.015 and fo = 10~ 84 (W m~ 2 ). This result implies 
that the magnitude of the flare (measured here with the GOES soft X-ray flux) 
is directly related to the free magnetic energy stored in the active region before 
the flare. If we include the X3.4 GOES flare of AR 10930 (Figure 8, case SI), 
we see that the observed free energy ratio qat = E^/Ep — 1 w 0.112 lies along 
the same trend of the extrapolated exponential function. 

Also the free energy Ef rcc , an absolute measure of the nonpotential energy, 
shows a correlation with the observed GOES flux, approximately following a 
powerlaw (Figure 8, bottom), 

//goes \ _ /-Efccc\ 3 / 17 a 
Wo J ~ V E, ) ■ [ > 
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Figure 8. Correlation between the free energy ratio gf rce = E{ IBe /Ep (top panel) or the 
free energy Ef Ice versus the logarithmic (preflare background-subtracted) GOES 0.5-4 A flux 
/goes for active regions A, B, C, and D, are shown, determined with the forward-fit method 
(diamonds) and with the misalignment method (crosses). Linear regression fits are indicated 
between linear or logarithmic values are indicated, yielding an exponential function (top panel) 
or a powcrlaw function (bottom panel) . 

with the constants / = 10 -8 4 (W m -2 ) and E = 6 x 10 30 erg. Only the active 
region (B) with the lowest free energy is an outlier to this powerlaw relationship. 

4. Discussion 

4.1. Accuracy of the Free Energy 

Compared with the other (numerical) NLFFF codes, our analytical NLFFF code 
has the following advantages: (1) computational speed that allows fast forward- 
fitting to observed coronal data; (2) simplicity of explicit analytical formulation; 
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(3) spherical geometry of solar surface is fully implemented; and (4) the free 
energy meets the five criteria of positivity, additivity of energy (or orthogonality 
of magnetic field components), positive scaling with force-free parameter a, 
potential field limit for small a's, and finiteness of nonpotential energy. In this 
study we calculated the free magnetic energy for four different test datasets and 
aim to validate the accuracy. 

The first test dataset consists of simulated loops, which are analytically de- 
fined by the same parameterization as the magnetic field lines that are forward- 
fitted to the simulated data. We find that the total nonpotential magnetic 
energies £"n are retrieved with an accuracy of « 10~ 5 for these cases, and the 
free energy ratios E-^/E-p are retrieved with an accuracy of m 10~ 3 . The high 
accuracy just confirms the convergence behavior of our code in the case of a 
unique solution, as it is the case when the forward-fitting model has the same 
parameterization as the fitted simulation dataset. 

The second test dataset consists of an analytical NLFFF solution, which has 
a completely different parameterization in terms of the Green's function applied 
to constant-a point sources (Chiu and Hilton, 1977; Low and Lou, 1990; Lothian 
and Browning, 1995; Malanushcnko et at, 2009, 2011, 2012). The particular case 
described in Low and Lou (1990) consists of three smooth elliptical magnetic field 
concentrations and contains an amount of 13% free energy, while we recovered 
only 2% with the forward-fit. The poor performance is possibly due to the partic- 
ular topology of the Low and Lou (1990) model, which has extended (elliptical) 
magnetic field regions parameterized with smooth Legendre polynomials, which 
cannot be fitted uniquely with spherical point sources. However, because this case 
is untypical for solar observations, which appear to consist of many small-scale 
flux concentrations that can be easier fitted with point sources, the accuracy of 
free energies may be much better for real solar data. 

The third data set consists of an observed flaring active region with a large 
X3.4 flare, for which a NLFFF solution was calculated. For the free energy we 
obtain a value of order 11%-18% before the flare, which approximately agrees 
with the most reliable NLFFF solution (16%-32%) computed in Schrijver et 
at, (2008) and Malanushcnko et at, (2012). Also for the sign of the change in 
nonpotential magnetic energy during this X3.4 flare we obtain the same sign and 
a similar value for the decrease of free energy (5% — 13%) as the values computed 
in Schrijver et at, (2008) and Malanushenko et at, (2012), « 13%. 

The fourth test dataset consists of four solar active regions, for which the 3D 
coordinates of coronal loops could be stereoscopically triangulated. For these four 
cases we obtained free energies of 0.6%, 2.3%, 4.4%, and 8.5%. Only for the first 
case (active region A, 2007 April 30) we have comparisons with other NLFFF 
codes. The most reliable NLFFF code that exhibits the smallest misalignment 
angle with stereoscopic loops, yield 3% free energy, which is close to our result 
of Efree/Ep = 0.6%. The free energy obtained from all 11 NLFFF solutions 
yields a much larger scatter, i.e., Ef ICC /Ep = 10% ± 12% (Table 1 in DeRosa et 
at, 2009). 
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4.2. Scaling Law between Magnetic Energy and Flare Soft X-ray Flux 

We found a correlation between the free magnetic energy and the GOES flux of 
the largest flare that occurred during the observing period (Equation 16). Since 
our nonpotential field solution is parameterized by a vertical twist of magnetic 
charges, the free energy is directly proportional to the magnetic energy associated 
with the azimuthal field component B v , and thus approximately proportional 
to the squared force-free a parameter or the number iVtwist of twists per length 
unit, (Equations (10) and (6)), 

feee(x) = -1=5 = (brsinef oc N 2 wist cc a 2 . (18) 

On the other side, the stress-induced heating rate En by Ohmic dissipation 
(or Joule dissipation) is proportional to the square of the currents, (e.g., van 
Ballegooijen, 1986), which is also proportional to the square of a (with j/47r = 
V x B = qB), and thus to the free energy, 

j 2 a 2 B 2 . . 

En = - oc , 19 

a Ana 

where a is the classical conductivity. Thus we can express Equation (16) in terms 
of twists, currents, or heating rates, 

Fsxn oc exp (q ircc ) oc exp (iV t 2 wist ) oc exp (a 2 ) oc exp (j 2 /B 2 ) oc exp (E^/B 2 ) . 

(20) 

Our empirical finding suggests that the energy radiated from the heated plamsa 
is not just proportional to the heating energy, but "exponentiated" by the current 
density j, the force- free parameter a, the number of twists A^wistj or the heating 
rate En of Joule dissipation. This implies a highly nonlinear mechanism that 
converts vertical twist into thermal radiation via dissipated currents, as envi- 
sioned in stress-induced reconnection (Sturrock and Uchida, 1981; Parker, 1983). 
Our linear regression fit of the free energy versus the GOES flux yields a lower 
limit of Fqoes = 10 -8 4 , which corresponds to a GOES-class Z4-event, about 
the magnitude of the smallest detectable nanoflarc. The largest flare among 
our analyzed cases is an X3.4 GOES-class flare. Apparently, the exponential 
relationship between soft X-ray flux and free energy discovered here holds even 
approximately for large GOES X-class flares (Figure 8). The nonlinearity is also 
reflected by the fact that the thermal energy of the soft X-ray flux is not just 
proportional to the free magnetic energy, but rather exhibits a highly nonlinear 
powerlaw function with a powerlaw index of « 3 (Equation (17)). 



5. Conclusions 

We calculated the total free magnetic energy contained in a coronal volume 
encompassing an active region for four different kinds of datasets: (1) simulated 
data, (2) data created from an analytical NLFFF solution by Low and Lou 
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(1990), (3) a flaring active region, and (4) four active regions observed with 
STEREO and SOHO/MDI. The free magnetic energy E frcc = E N - E P is de- 
fined by the difference of the nonpotcntial (-En) and the potential magnetic field 
energy (-Ep). The nonpotcntial magnetic field Bn(x) is defined by an analytical 
approximation of a NLFFF solution that is parameterized by buried magnetic 
charges with vertical twists (derived in Paper I). The numerical code that per- 
forms fast forward-fitting of magnetic field lines to coronal 3D constraints, such 
as stereoscopically triangulated loops, is described in Paper II, along with the 
simulated data. Our findings are as follows: 

1. A first method to calculate the free energy results from forward- fitting of our 
analytical NLFFF approximation by associating the perpendicular magnetic 
field component B± with the free energy, dEf Iee = B"j_/8w, while the parallel 
component B\\ is associated with the potential field energy dEp = B^/8ir. 
This definition of the free energy fulfills the conditions of (i) positivity of free 
energy, (ii) additivity of energies, E^ = Ep+E[ rce and orthogonality B±(x) _L 
£>p(x), (iii) a positive scaling with the force- free parameter, dE{ 100 oc a 2 , (iv) 
the potential- field limit, E-^(a M> 0) = Ep, and (v) the finiteness of the 
nonpotcntial energy with height h, E^(h n> oo) = fi max . 

2. A second method to estimate the free energy can be obtained from the mean 
misalignment angle A/i = /ip — fx^ between a potential and a nonpotential 
field (or stereoscopically triangulated coronal loops). The free energy ratio 
is then E{ ree /Ep m tan 2 (A/x). Wc find that the uncertainty of this method 
amounts to « ±2% for the nonpotcntial magnetic energy. 

3. Calculating the free energies for the simulated data we find a high fidelity 
of order 10 -5 in retrieving the free energy, which is due to the fact that the 
simulated data have the same parameterization as the forward-fitting method, 
constraining a single best-fit solution. 

4. Calculating the free energy for the Low and Lou (1990) analytical case, our 
NLFFF code finds a significantly lower value than theoretically calculated, 
probably because of the special morphopology (parameterized with smooth 
Legendre functions), which cannot adequately be fitted with our NLFFF 
code that is designed for spherical magnetic point sources, as found in solar 
magnetograms. 

5. Calculating the free energy for observed active regions constrained by the 3D 
coordinates of stereoscopically triangulated coronal loops, we find free energy 
ratios of qf ICC = Ef Tee /Ep rts 1% — 10%. The uncertainty of the free energy 
determined with our forward-fitted NLFFF approximation appears to be at 
least as good as the uncertainty among other (standard NLFFF extrapolation) 
codes. 

6. We find also a correlation between the free magnetic energy -Ef roo and the 
GOES flux of the largest flare that occurred during the observing period, 
which can be quantified by an exponential relationship, Fgoes oc exp(<7f ree ), 
implying an exponentiation of the dissipated currents. 

In summary, this study demonstrates that the free energy in active regions can 
be calculated and predicted with our analytical NLFFF approximation with an 
accuracy that is commensurable with other standard NLFFF codes. Our code has 
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the additional advantages of computational speed for forward-fitting of coronal 
data, correct treatment of the curved solar surface, positivity, and finiteness of 
free energy. In addition, forward-fitting of our NLFFF approximation achieves 
a significantly smaller misalignment angle with respect to the observed coronal 
loops (fj, rj 2° — 22°), compared with the results of other NLFFF codes (fj, = 
20° — 44°; DeRosa et al, 2009). The most limiting drawback of our method is the 
availability of stereoscopic data with suitable spacecraft separation angle (which 
was most favorable in 2007, the first year of the STEREO mission). In future 
work we attempt to circumvent the 3D geometry of coronal loops by using only 
the 2D projections of coronal loops, which can (manually or automatically) be 
traced from loop-rich EUV or soft X-ray images and do not require stereoscopic 
data at all. 
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